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The origin of “secular” behavior for the nonlinear cold electron plasma equations is studied. The X V 

equations involved are closely related to the Klein-Gordon equation with a small nonlinear term. 

A method is developed for arriving at perturbation theoretic solutions of this equation, and the 
method is then applied to the case of the higher-order effects of an electromagnetic wave propagating 
in the cold electron plasma. An explicit expression for the second order frequency shift is calcuIatecL^^^^^^ 


n 




I. INTRODUCTION 

HEN one attempts to calculate the behavior 
of nonlinear mechanical systems in pertur- 
bation theory which goes beyond the linearized ap- 
proximation, often (though by no means always) 
the difficulty of “secular” behavior appears. The 
higher orders contain, in addition to trigonometric 
terms, time-proportional, or “secular” terms. The 
unbounded character of these terms soon invalidates 
the perturbation theoretical assumptions of smallness 
on which they were derived. 

A systematic program for doing a type of perturba- 
tion theory which is free of secular terms in the case 
of the harmonic oscillator equation with a “small” 
nonlinear term was given some time ago by Krylov 
and Bogoliubov, 1 and later refined and mathe- 
matically justified by Bogoliubov and Mitropolskii. 2 
Recently, considerable interest has arisen in modi- 
fying these techniques to deal with partial dif- 
ferential and differentio-integral equations, partic- 
ularly in connection with the work of Frieman 3 and 
Sandri. 4 Earlier calculations, which are more closely 
related to this work, were made by Jackson 5 and 
Sturrock. 6 

The system treated below is a partial differential 


* Permanent Address: Institute for Fluid Dynamics and 
Applied Mathematics, University of Maryland, College Park, 
Maryland. 

1 N. Krylov and N. Bogoliubov, Introduction to Nonlinear 
Mechanics , translated by S. Lefshetz (Princeton University 
Press, Princeton, New Jersey, 1947). 

2 N. Bogoliubov and Y. A. Mitropolskii, Asymptotic 
Methods in the Theory of Nonlinear Oscillations (Gordon and 
Breach Science Publishers, New York, 1961) [translated from 
the Russian]. A quite similar technique was put forward 
independently by M. J. Lighthill [Phil. Mag. 40, i 179 (1949)]. 

3 E. A. Frieman, J. Math. Phys. 4, 410 (1963); and (to 
be published). 

4 G. Sandri, Ann. Phys. (N. Y.) 24, 332 (1963). 

6 E. A. Jackson, Phys. Fluids 3, 831 (1960). 

6 P. A. Sturrock, Proc. Roy. Soc. (London) A242, 277 
(1957). See also: C. L. Dolph, J. Math. Anal. & Appl. 5, 
94 (1962); H. S. Tsien, in Advances in Applied Mechanics 
(Academic Press Inc., New York, 1956), Yol. IV, p. 281. 


equation, considerably simpler than those described 
in Refs. 3 and 4, but which we nonetheless believe 
to be of value in illuminating, in a relatively un- 
cluttered way, some of the new features which 
emerge as a consequence of the fact that we are 
studying a partial differential, rather than an ordi- 
nary differential system. 

The equation treated by Bogoliubov and Mi- 
tropolskii is 

(& + 4«> - 4' I) • <•> 

where e is a formal expansion parameter (eventually 
to be set equal to one) used to indicate the relative 
“smallness” of the right-hand side. F(x, dx/dt) is 
a known nonlinear functional of x and dx/dt ; co 2 
is a constant. 

We examine the equation 

where e is again the formal expansion parameter, 
F is again a known functional (which depends on 
the problem under consideration), and c 2 and A 2 
are non-negative real constants. 

The immediate motivation for studying (2) was 
that the equations of motion for a “cold” electron 
plasma in a uniform positive background may be 
put in this form, component by component. How- 
ever, the cold plasma equations are superficially 
less simple than Eq. (2), so we shall consider it 
first. 

We shall first give an example of both secular 
and nonsecular behavior for Eq. (2), and show how 
the difficulties may be remedied in the secular case. 
Then in Sec. Ill, we shall give a simple physical 
example, that of the second order behavior of a 
nonlinear electromagnetic wave in a “cold” elec- 
tron plasma. 
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COLD PLASMA EQUATIONS 


II. THE KLEIN-GORDON EQUATION WITH A SMALL 
NONLINEAR TERM 


A. Single Monochromatic Wave 

We first consider perturbations about the follow- 
ing 6 = 0 solution to (2) 


/ = a cos (K 0 x — o)J + <t>) f 

(3) 

where 


2 2 7^2 I ^ 2 

w 0 = c K 0 + X 

(4) 

and for 6 = 0, the quantities K 0 , co,„ a, 

and <f> are 


constants. 

This happens to be a situation for which secularly 
arises, so we shall set up the necessary formalism for 
handling it from the beginning. Following Ref. 2, 
we seek a solution to Fq. (2) of the form 

/ = a cos \p + eu^a, i) + t 2 u 2 (a, \f>) + • • • , (5) 

where the amplitude a is now determined as a 
“slowly varying” function of x and t by the relations 

da/dt = eiti(a) T - e 2 A 2 (a) T - • ■ • , (ha) 

da/dx = e Di(a) + e 2 D 2 (a) + • • • . (Gb) 

^ is a new “phase” variable, to be chosen to coincide 
with the phase of (3) for e = 0, 

d\[//dt = — co 0 + eB l (a) + e 2 B 2 (a) + • • • , (7 a) 

dt/dx = K 0 + eC\(a) + e 2 C\(a) + ■ • • . (7b) 

The (as yet undetermined) functions A u A 2 , • ■ ■ , 
Ri, B 2 , * • • , C\, C 2 , ■ • • , D lf D>, ■ • • , are to be 
chosen so as to render the solution (5) free from 
secular terms. Here, “secular” must be interpreted 
to mean ^-proportional, i.e., solutions can break 
down due to linear growth in x as well as t. The 
functions u X} u 2 , • • • , are to be periodic in \p. 

The program is to express the various terms in (2), 
by means of the relations (5), (6), and (7), in terms 
of the u’ s, A’ s, B’s, C’ s, and D 1 s, as functions of 
a and \p. For brevity's sake, we shall first only go to 
0(e), though in principle the method may be carried 
to any order in 6. Generally, the algebra becomes 
prohibitive beyond 0(e 2 ). Hereafter, the notation 
“+ • ♦ •” will mean “of higher order in e.” 

The result of differentiating / with respect to x 
and t and using (6) and (7) is 

^rL — c 2 ^4 + X 2 / = a cos " <x>l + c 2 Kl + X 2 ) 
at ox 

+ + «.) + 2 KA, + <?K 0 D t ) sin f 


The zeroth order part of (8) vanishes identically, 
by virtue of (4) . The coefficient of e is to be equated 
to F(fj df/dt, d//d;r), with the arguments replaced 
by their zeroth order values: 

/— >a cos 4 / , d//cM— >co 0 asin \p, d//d£— > —K 0 asm yf/. 

For the purpose of finding u lf it is most convenient 
to write F as a Fourier series in 4 / 

F(a cos 4 b sin 4 / , — K Q a sin 4 /) 

00 

= 0o (a) + [0„(a) cosnt/' + /„(a)sin»£|. (9) 

n= 1 

The f n (a ) and g n (a) are known functions of the ampli- 
tude a , determined by the functional form of F. 
We shall make the assumption, always satisfied in 
practice, that the f n (a) and g n (a) go to zero, as 
a goes to zero, at least as fast as 0(a). 

The equation for u x may be written as 

\\d 2 /dt 2 + IK = g 0 (a) 

+ [A (a) - 2u ( >A , - 2c 2 K 0 D x ] sin ^ 

+ [<7i(a) — 2 au {) B x — 2 ac 2 K 0 C x ] cos ^ 

+ 12 \Qn(a) cos 71 4* + /»(a) sin n\p ]. (10) 

n-2 

This is, effectively, an ordinary differential equa- 
tion in 4* of a standard type. Its solution contains 
terms proportional to 4 / sin 4 / and ^ cos 4 / J unless 
the coefficients of sin 4 / and cos 4 / on the right hand 
side of (10) vanish. Jj they vanish, then one may 
find the general solution in the form 

u x = v {) (a) + 12 [v n (a) cos + w„(a)s mnyp]- (11) 

Our obj (active is to find a solution free from \k- 
proportional terms. This impels us to choose 
B i, Ci, and D x , so that 

2(« 0 il 1 + c 2 K 0 D x ) = A(a), (12a) 

2a(o ) 0 B X + c 2 K {) C x ) = gi(a). (12b) 

Two more relations may be deduced from the 
conditions 

d 2 4'/dt dx = d 2 i^/d.r dt , d 2 a/dt dx ~ d 2 a/dx dt. 

Reference to Eqs. (6) and (7) shows that, to lowest 
significant order, these become 

A x dCi/da = D l dBJda, (13a) 

A x dD x /da = D x dAi/da. (13b) 

Equations (12) and (13) are four relations for 
the four unknowns A Xj B u C i, and D if and once 
have been obtained, (10) can be solved in a 
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straightforward way by means of (11), and u x will 
contain no ^-proportional terms. can be inter- 
preted as a “frequency shift” and C x as a “wave 
number shift.” 

However, it is impossible to solve (12) and (13) 
completely without specifying the physical problem 
in more detail. For instance, if we wish to work a 
boundary value problem in which f is required to 
oscillate sinusoidally at a given x for all t, we do 
not expect a or 4 / to vary with t , and B x and A x 
may be set equal to zero. This leaves 


Cl = g l {a)/2cK u a 

(14a) 

as the “wave number shift” and 


/>, = «a)/2cX. 

(14b) 


If we are interested in an initial value problem 
in which a pure sine wave is given for all x at t ~ 0, 
we may correspondingly set C x and D x equal to 
zero and solve for A x and B x . Other choices are 
necessitated by still other problems. (Note in passing 
that in both cases, the corrections — > 0 as a — > 0, 
as they must in order to make physical sense.) 

All this has been for a completely general F, 
imagined to contain ail harmonics in 4 / - In practice, 
it often happens that F contains only a few har- 
monics. Observe that secularity may not arise for 
some forms of the nonlinear term. Thus, if F is 
proportional to (/) 2 or / df/dx, say, f x (a) and g x (a) 
will both vanish, and in this order, we may set 
A i, B i, Ci, and D x , all identically zero, which is 
equivalent to a completely straightforward kind of 
perturbation theory. On the other hand, if we were 
to have F proportional to (/) 3 , f x (a) is still zero, 
but <7, (a) = 3a s /4 and a straightforward perturba- 
tion theory no longer works. The occurence of 
secularity thus depends, in a given order, on the 
particular form of the nonlinear term. 

To close this discussion, we give the full solution 
for u x (a, \p) of Eq. (11). If we make the (arbitrary) 
choice that all the first harmonic shall be collected 
in lowest order, we may choose v x (a) — w x (a) = 0. 
Then for n ^ 1, • 

*>»(«) = gn{a)/\ 2 {-n~ + 1 ), ^ 15a ^ 

«>»(«) = /„(a)/X 2 (-n 2 + 1), 


'fn(a) ! 

-If 

2 IT 

dyp 

sin n\p 

g«(a) * 



cos n4 e 


•F(a cos 4 >, co 0 a sin \p, —K 0 a sin ^). (15b) 


When the phenomenon of secularity occurs, one 
must go to one higher order in e in order to get 
corrections to / which are uniformly valid through 
terms of 0(e) for changes in x and t of 0(l/e), and 
so the solution (15) is not of great interest by itself. 2 
However, the expressions such as (14) for frequency 
and wave number shifts are accurate £o 0(e), and 
are usually more accessible to measurement than the 
expression for u Xl in any case* 

B. Two Coupled Monochromatic Waves of Differing 
Frequencies 

AVe now proceed to the case in which (2) has the 
e = 0 solution 

/ = a[cos (K x x — u x t + <f> x ) 

+ cos (K 2 x — o o 2 t + <l > 2 )] , (16) 

where 

co? - c 2 K\ - X 2 = co? - c 2 K\ - X 2 = 0, (17) 

but where there is not necessarily any other relation- 
ship between coi, K x and co 2 , K 2 . If we were to do a 
straightforward perturbation theory, assuming 

/ = / (0) + e/ (I) + e 2 f V2) + ••• , (18) 

with / (0) given by (16), it is clear that we can always 
find a secularity-free perturbation theoretical solu- 
tion to (2) of the form 

= Y /'"V) exp [i[r(KiX — to , t + <f>,) 

r ,a = — co 

+ s(K 2 x — C 0 2 t + <£ 2 )]} (19) 

except when there happen to exist two integers r, $ 
such that 

X 2 (-f 2 - s 2 + 1) - 2rS(anco 2 - c 2 IC x K 2 ) = 0. (20) 

When (20) is fulfilled, / (1) , like the u x of Eq. (10), 
will in general contain terms proportional to the 
phases, as well as trigonometric terms. 

For a randomly selected pair of values to,, K ,, 
and co 2 , K 2 , there will exist no integers ?, s which ful- 
fill (20), and there is no necessity for the Bogoliubov 
techniques. When, on the other hand, it happens 
that two such integers r, s, do exist [the relation (20) 
will also hold for — f, — s, of course], then the 
Bogoliubov techniques are called for. One method 
of introducing them is the following, though it is 
not the only method. [ Note added in proof. For a 
completely general F(f } df/dt , df/dx), there will 
exist nonvanishing Fourier coefficients on the right- 
hand side of the equation for / (n) corresponding to 
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terms in which one of the integers r, s is zero and 
the other is plus or minus one. For this case, Eq. 
(20) is, of course, trivially satisfied, and the phe- 
nomenon is only a slight generalization of that of 
the previous subsection. For the point made here, 
we shall assume that F is of such a nature that, to 
the order we are going, there exists only one pair 
of integers r, s satisfying (20), neither of which is 
zero, and for which the Fourier coefficient corre- 
sponding to these integers on the right-hand side 
is nonzero. We are indebted to Dr. J. Iv. Hale for 
a conversation which led to the discovery of an 
error on this point in our original manuscript.] 
Seek a solution of the form 

/ = a (cos \p ! + cos ^ 2 ) + m x (a, $ u \p 2 ) + * • * , 

where now the amplitude a and the phases x and ^ 2 
develop according to 

da/dt = eA(a , ^ 1 ? ^ 2 ) + * * * , 

da/dx = eB(a , \f/ 2 ) + * * * , 

d\p u2 /dt = —co lr2 + eC li2 (a, \f/ u $ 2 ) + ••• , 
^^ 1 , 2 /d.T = K lt2 + dJi i2 (a, ^ 1 , ^ 2 ) + *•* . (21) 


The zeroth order part of (22) vanishes identi- 
cally. The first order part now must be equated to 
eF(a( cos + cos yp 2 ), a (an sin + a> 2 sin \p 2 )j 
— a(Ki sin i/'i + K 2 sin ^ 2 )). It is most convenient 
to write this as a complex Fourier series, 

F = E (23) 


where F* n = since we deal only with real 

quantities. Calling everything in the first curly 
bracket in Eq. (22), (7(a; .4; B\ C i (2 ; D 1(2 ), we may 
finally write the equation for u A as 


x2 (fi + fi + Ui ) + 2(w “° 2 ~ c2RiK2) 

= G(a ) A ) B ,* Ci >2 \ D \ >2 ) 

+ 00 

+ E (24) 


m , 7i = — a> 


Equation (24), analogously to Eq. (10), has a 
readily obtainable, secularity free solution of the 
form 


The left-hand side of (2) may again be computed, 
using the relations (21), to give 

— — c + X 2 / = u[( ci) x + c 2 K'i + X 2 ) cos \p x 
ox 


dt 


4" ( — w 2 + c 2 K 2 -f- \“) cos ypz\ 

I / . , To . I dC, , dC 1 

+ sin 2o)iA + dUi + ao> 2 

l L d^ x d\p 2 




+ 2 c 2 K,B + ac 2 K, ^ + ac 2 K dDl 
d\p x 

+ sin i p 


d^ 2 J 

O A I I dCo 

Z(jl 2 A -f- (ZCOx ~ ■ ' GC0 2 


dC 2 

dipi 


dC, 

d\[/ 2 


+ 2c 2 K,B + acK, + ac 2 K, 


dti 


+ COS \p 1 


ai> 2 

d\p 2 

dA 


2au i C l — Ui ~~ — co 2 ■ . 

(Jyi d\f/ 2 


dB 


n 77- .2 7/- d.B 


+ 2ac K 1 D l - c 2 K, - c 2 IC 
d'Pi 


+ COS i/": 


o dA 

Zdo) 2 O 2 w x w 2 - . 

0^2 


aB 


+ 2c 2 aK 2 D 2 - FK, yj- - c 2 K 2 

dy 1 d^ 2 J 

, j.Je 2 u, d 2 Ul \ 

+ V v#! + ^! + N 


+ 2(«i&) 8 — c 2 KA\ 2 ) 


d\ \ 
dfi dtJ 


+ 0 (< 2 ). ( 22 ) 


uM u a) = E (25) 

to ,n- — » 

if and only if the functions A, B } C\, 2 , D l 2 are 
chosen to depend on and \p 2 in such a manner 
that: 

^ d\p- 2 exp [—i(fip i + S^)] 

•G(a; A ; B; C,.,; D llt ) - ^,.-( 0 ) = 0. (26) 

Since (26) involves complex numbers, it really 
amounts to two real equations upon A } B } Ci, 2j 
D 1(2 , which are linear partial differential equations 
with periodic coefficients. (We may consider A, B, • * • 
as being expressed as Fourier series in yp u ^ 2 ; in 
this case, (26) leads to algebraic equations for the 
Fourier coefficients.) We need not write them down 
in full detail, since they are most cumbersome, and 
not of interest for our purposes here. 

Three additional conditions may be adduced from 
the requirements that d 2 a/dx dt — d 2 a/dt dz, 
dVi, 2 /d.r dt = dVi, 2 /d£ dx; they are 
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Satisfaction of these five relations on the six 
quantities A, B , C lf2 , D u2 , then, will guarantee 
that the solution shall be secularity free. It is clear 
that considerable latitude is left to choose the func- 
tions conveniently for whatever problem is under 
consideration. Nor does secularity necessarily have 
a simple physical interpretation as it does for Eq. (1), 
where the presence of secular terms is simply inter- 
pretable as a resonance between one of the fre- 
quencies present in the nonlinear coupling term and 
and one of the natural frequencies of the system. 
It would take considerably more insight to appre- 
hend the physical meaning of Eq. (20), for the case 
in which neither f nor s is zero, or any connection it 
might have with the secularity condition for Eq. (10). 

III. THE COLD ELECTRON PLASMA 

Suppose we consider a cold (no thermal motions) 
electron plasma of equilibrium number density n 0} 
moving in a uniform positive background, assumed 
immobile. If we make a perturbation about a uni- 
form, field-free equilibrium, the appropriate vari- 
ables for describing the system are v, the electron 
velocity, — e[n 0 + n(x, Z)L the electron charge 
density, +en 0 , the positive background charge den- 
sity, and the electric and magnetic fields E, B. All 
these will be* treated as perturbations — i.e., first 
order in the amplitude — except for n 0 . 

The dynamical equations are well known: 



dn 

dt + 

+ 

© 

II 

V 

o 

(28a) 

dv 

dt 

+ V*“" 
dx 

m N 

E + Jvxb) , 

(28b) 



1 dB _ 

c dt 

j-xE, 
dx ’ 

(28c) 


idE _ 
c dt 

= S* B + 

(no + n)v; 

(28d) 



II 

W 

*1* 

-4tt en, 

(29a) 



dx 

= 0. 

(29b) 

Equations (29) can be regarded as initial conditions; 


once fulfilled, they are preserved by Eqs. (28). 

A small amount of algebraic juggling shows that 
Eqs. (28) lead to 


f?* + c ' - -<*« 

(""(’■£ v+ = txB ) • 
d ( d dB\ , : SB 

h (OT >] ■ 


d 2 dB , 2 

dt 2 at + c 


— 47 rec — x - 
dx 


(30a) 


(30b) 


where u>l e = 47r n 0 e 2 /m is the plasma frequency. The 
formal expansion parameter e has been written on 
the right-hand side of Eqs. (30) only to remind us 
that these terms are “small” in the sense of being 
second order in the amplitude. 

If we now restrict ourselves to disturbances which 
are functions of only one spatial dimension (x, say), 
we may write the expression (30a) in the form 

(!?- ^ I? +^+ < 3i > 


where 5 ui = 0 unless i = 1 and 5 1(1 = 1. gq repre- 
sents the fth component of the right-hand side of 
(30a). Since Eq. (30b) is just the curl of (30a), it 
is not necessary to write it in the form of (31). 
It matters not a bit that we have not bothered to 
express the right-hand side of (30a) as a function 
of E alone, since we always need the values of the 
quantities which enter into to one lower order 
than those under consideration on the left-hand 
side of (30a), and these will be shown to be ob- 
tainable in each order directly from (28). 

The 6 = 0 solution to (31) is just the standard 
cold-plasma set of field-free normal modes, which 
is well understood. 7 It will be apparent that each 
component of (31) has the form of Eq. (2) and is 
therefore immediately susceptible to the methods 
of Sec. II. Equation (31) represents a generalization 
of the equation treated by Jackson, 5 which, however, 
does not exhibit all the features of (31), due to the 
absence of ^-derivatives on the left-hand side in 
the purely electrostatic case, and the fact that 
the x and t dependences separate. 

The number of possibilities from (31) is very 
large, due to the wide range of choices for the e = 0 
solution. The simplest case (beyond that of a pure 
electrostatic oscillation 5 ) is that of a pure transverse, 
linearly polarized, electromagnetic wave. We shall 

7 I. B. Bernstein and S. K. Trehan, Nucl. Fusion 1, 3 
(1960). 
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extend this solution to the next two orders above 
the linear approximation in the following subsection. 


A. Plane Electromagnetic Wave 

Calling our basis vectors e x , $ V) e z , the appropriate 
zeroth order solution is 


E = E 0 e v cos ( k 0 x — o> Q t), 

B = — EA cos (ka - o) 0 0, (32) 

cE 

v = — e u sin (k 0 x — co o 0, n — 0, 
rm o 

where 

cog = a, 2 e + C 2 kl (33) 


This zeroth order solution for E is conveniently 
written: 

E t = + <f>%e~ iikoX " Uot) ] 

where $ is a column vector 

fo] 


loj 

and a — E 0l the amplitude. 

The results of Sec. II suggest that we seek a 
solution of the form 


E t = a(4>,e" ! ' + *V") + eitf'ip, #) 

+ e 2 «; 2) (a, \p) + • • • , (35) 


where now the amplitude a and phase ^ are to vary 
in x and t according to 

da/dt — cA j (a) 4“ e A 2 {ci) -b • * • , 

da/dx = eDi(a) + e 2 D 2 (a i) + • ■ ■ , 

dip/dt — — w 0 "b cJ5i(a) + e 2 B 2 (a) + • * * , ^ ^ 

d\p/dx — k 0 + €(7i(a) + eC 2 (a) + • * * , 


with the functions A l} B 1} • • * , as yet undetermined. 

The left-hand side of (31), expressed in terms of 
these functions, becomes: 
for i — 1, 


(f? + 03 °°) Ei ~ 4 


2 O Xl\ 
co 0 


dr 


+ &pe U l + 2( — io2 Q A x 


~b cucoBi)^^' * + 2(ia> 0 ^4i + ao>oBi)<t>fe 


*0-' + 



+ 2(— iu 0 Ai + aw i) B l — ic 2 k 0 D 1 + ac 2 C 1 k 0 )<f) l e" 1 ' 

~h 2(io) 0 A l + ao) 0 Bi + fc 2 /c 0 Di (38) 

+ + o(/). 

These three expressions are to be equated to the^ 
three components of the right hand sides of Eqf. 
(31), with the quantities v, E, B, n 2 'eplaced by tl?eir 
values (32), and with the substitution E 0 -*> a. 
The result can be conveniently written as a Fourier 
series 


5, = £ *,(ro, a)e im \ (39) 

ffl** — CD 

\ 

where the <f\(ra, a) are known functionals of a. 

Suppose we now consider whether secular terms 
can arise. We seek the solution for u\ l) in the form 

a)e tm * . (40) 

— co 

The equation for u\ 1] then becomes 
£ [(-m 2 + l)c4 + m 2 S 1 . i (a o 2 pe - co 2 )]t u< v (m, a)e’“* 

til =* — 03 

==: 2( zcoqAj - b clo)qB i 'ic JcqDi ~b uc C\ko)(j)i6 

— 2(iu Q A x + acooi?! + ic 2 k 0 D l + ac 2 C\k () )<j>*e~ l4 
+ (41) 

tn *= — as 

(we have made use of 4>\ — 0). 

We may always solve (41) for u™(m, a) with the 
absence of ^-proportional terms, except when a 
coefficient of i4 l) (ra, a) on the left hand side of (41) 
vanishes for some m. Since o % > co 2 e it is clear that 
this can only occur when m = dbl, and only in 
the components i = 2 or 3. However, a term-by- 
term inspection of <f\(ra, a ) reveals that for i ~ 2 
or 3, the quantities $i(m, a) vanish identically for 
the choice (33) of the zeroth order values. Thus a 
result emerges which would have been hard to guess 
from the original set of equations (28): secular be- 
havior cannot arise in second order for this case, 
and we may set A x = B x = C x = D x = 0, and 
solve in a completely straightforward way for 
t4 n (m, a) 

u\ v (m, a ) = 4> t (m, a)[(— ?n 2 + l)a>p e 

+ ~ o>l)T\ 

a) = 0. 


for i — 2 or 3, 


+ 0(e 2 ); 

(37) 


m 9^ dhl, 
(42) 
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The calculation of the coefficients of the Fourier 
coefficients is a matter of simple algebra, and we 
may write down in full the first correction to the 
linear solution (32) : 


u 


(i) 


n 


T-'d) t,i 2 cfcpCOpe sin 2(k 0 x ccpt) 

* ~ 0 2nu4 (3c0p e + 4 c 2 kl) 

B (1) = 0, 

e 2 n 0 kl cos 2 (k 0 x — co 0 t) 

- 2 - 2 (3(4. + 4 c 2 kl) ’ 


v (l) = 


-m 2 2 

m co 0 

y T 2 6 fcg COS 2(fco^r ^0^) a 

m co 0 (3to pe + 4c /c 0 ) 


The only qualitatively new feature which shows 
up in 0(e), then, is a longitudinal electric field. It 
is only in 0(e 2 ) that secularity manifests itself and 
it becomes necessary to use the Bogoliubov methods. 
Making use of the 0(e) solution we have just com- 
puted, we may write down the 0(e) part of the 
left-hand side of (31) in terms of the A’ s, B’s, etc. 
We need 


I d 2 E { 
\ dt 2 


a d 2 E { 


+ co p JSi i 

OX ) O (t*) part 


32 ( 2 ) 


2J / 2 2/2\ d Uj 2 <2) 

= e y<j)Q — c k 0 ) — 2 ~ + 0)pe^» 

~T ( — 2fc0o'd2 + 2q<j$qB 2 2ik 0 c T ) 2 


T 2ac 2 /c 0 C 2 )</> t -e ^ (2iajoA2 2aojQB 2 
+ 2 ik 0 c 2 D 2 + 2ac 2 fc„C' 2 )0*e~ i 4- 


(44) 


It is also a simple matter to compute the 0(e) 
part of the right-hand side of (31) from (43) and (32). 
The result is, expressed in terms of a and 


2 7,2,2 3 

^ -1 tO pe'^O*' U Qt, 2 

^ Jo(<3)part “ 4co 2 w 2 (3co 2 e + 4c 2 fc 2 ) 

• [3(e 3t * + c" 3 ^) - (e** + e'**)]. (45) 

Since the only nonvanishing component of (45) 
is in the 2-direction (y-direction), the only non- 
vanishing component of u\ 2) will also be u { 2 2) . The 
equation to be solved in the 0(e) approximation is, 
therefore: 


aV 2) 

(<4 - c 2 kl) + c olu™ = -(-2i» 0 A 2 

-f* 2atcoB 2 — 2ik$c E) 2 + 2afc 0 c C/ 2 )<f> 2 ^ 

— (2fw 0 A. 2 + 2aoooB 2 + 2ik 0 cD 2 


+ 2ac 2 k 0 C 2 )<t>%e‘ i ' t ' + 


2t223 
Wpe/Co^ & 


+ icX) 

■ [3(e s ‘* + e- 3,> ) - («** + «“'*)]. 


We may find a solution to (46) of the form 
u?\a, *) = £ ^ 2> (a, rn)e im \ (47) 


which is free of ^-proportional terms, if and only 
if the coefficients of e x * and e~ x * on the right-hand 
side of (46) vanish. This gives us two conditions on 
A 2 , B 2 , C 2} D 2 , which are algebraic relations. Two 
more, analogous to Eqs. (13), are given by the 
requirements that d 2 a/dt dx = d 2 a/ dx dt and 
d 2 yp/dt dx — d 2 \p/ dx dt; they are 

D 2 (dA 2 /da) = A 2 (dD 2 /da ), (48a) 

D 2 (dB 2 /da) = A 2 (dC 2 /da). (48b) 


Let us now specialize the problem to one in 
which the spatial periodicity and amplitude are 
given, and we are to calculate the “frequency 
shift.” This means setting C 2 = D 2 = 0. The con- 
dition that the coefficients of e 1 *, e~ x * on the right- 
hand side of (46) vanish reduces to the equation 


( ioiftAz * cuoqB^i) 


o>lXe 2 a 3 

4m 2 co 2 (3co 2 e + 4c 2 /c 2 ) 


(49) 


and its complex conjugate relation. We may there- 
fore find A 2 = 0, and the frequency shift Aco be- 
comes (substituting E 0 for the amplitude a) 


Aco = — B 2 


eWpul El 
4m 2 c4 (3a) 2 e + 4c 2 * 2 ) ’ 


(50) 


It does not seem worthwhile to write down the 
0(e 2 ) corrections to E, B, v, and n, though it would 
be easy to do so. In any experiment one might 
imagine, Eq. (50) would probably be the easiest 
quantity to measure. 


B. Standing Wave* 

It is possible to calculate the frequency shift of 
a standing electromagnetic wave of given periodicity 
as well. One assumes a linearly polarized standing 
wave solution for (31) in lowest order and writes 

Ei = cq sin k 0 x cos \j/ 

+ eu ( i l) (a ) ypj x) + eV 2) (a, x) + 

Then, only the ^-dependence of the phase variable 
and amplitude are expanded: 

da/dt = eAda) + e 2 A 2 (a) + • * ■ , 

d\p/dt — o) 0 + eB^a) T* e 2 B 2 (ai) 

In lowest significant order, the frequency shift turns 


(46) 


This section added in proof. 
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out to be exactly one eighth of the result given in 
Equation (50). 

IV. DISCUSSION 

We have given a technique for obtaining uni- 
formly valid, perturbation theoretic solutions to the 
Klein-Gordon equation with a small nonlinear term. 
The method has been applied to calculate the second 
order frequency shift of an electromagnetic wave in 
a cold electron plasma. The smallness of the expres- 
sion (50) for attainable parameters is an indication 
of just how good an approximation the linear theory 
is at these frequencies. 

It is not to be inferred, however, that the method 
adapts itself readily to all partial differential equa- 


tions. For instance, the reader can easily convince 
himself that it fails for a nonlinear sound wave de- 
scribed by the Euler equations. [It appears to fail 
in all situations for which the e — 0 equation is a 
wave equation, ( d 2 /dt 2 — cd 2 /dx 2 )f = 0.] The 
physical reason is that, due to a steepening of the 
exact nonlinear wave front — obtainable from the 
Riemann invariants — a vertical tangent develops 
after a time of 0 (1/amplitude). This destroys any 
regularity properties which may have existed in the 
original wave profile. Since the exact solution does 
not remain “close” to the e = 0 solution in any 
sense, after a time of order 1/6, it is not surprising 
that perturbation theory is of little use beyond 
this time. 


